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1 INTRODUCTION 



1 Introduction 

The problem of Maxflow was formulated by T.E. Harris in 1954 while studying the Soviet 
Union's railway network, under a military research program financed by RAND, Research 
and Development corporation. The research remain classified until 1999. The Maxflow 
problem is defined on a network which is a directed graph together with a real positive 
capacity function defined on the set of edges of the graph and two vertices s, t called the 
source and the sink. A flow is another function of this type that respects capacity constraints 
and a Kirchoff 's law type restriction on each vertex except the source and sink. The net flow 
of a flow is defined as the amount of flow leaving the source. The problem of maxflow is to 
find a fiow with maximum net fiow on a given network. In the first section we will define 
clearly such concepts and present basic results in the subject. 

Throughout the second section, we will present three different algorithms for the solution of 
Maxfiow. In 1956 L. Ford and D. Fulkcrson devised the first known algorithm that solves the 
problem in polynomial time. The algorithm works starting with the zero fiow and finding 
paths from source to sink where flow can be augmented preserving the flow and capacity 
restrictions. We then analyze a more efficient algorithm developed by A. Goldberg and E. 
Tarjan in 1988. This algorithm works in a different fashion starting with a preflow, a function 
saturating edges adjacent to the source, and then pushing excess of ffow to vertices estimated 
to be closer to the sink. At the end of the algorithm the preffow becomes a ffow and in fact, 
a maximum flow. Finally we describe Dorit Hochbaum's pseudoflow algorithm, which is the 
most efficient algorithm known to day that solves the Maxflow problem. 

In the third section we show the usefulness of this subject and present three applications of 
the theory of Network flow. First we show how well-known theorems in combinatorics such as 
the Hall's Marriage theorem can be proven using Maxflow results. We then show how to flnd 
a set of maximal chains in a poset with certain properties using the results in the previous 
sections. Finally we describe an algorithm for image segmentation, an important subject in 
computer vision, that relies on the relation between a maximum flow and a minimum cut. 

The problem of Maxflow is a widely developed subject in modern mathematics. Efficient 
algorithms exist to solve this problem, that is why a good generalization may permit these 
algorithms to be understood as a particular instance of solutions in a wider class of problems. 
In the last section we suggest a generalization in the context of simplicial complexes, that 
reduces to the problem of Maxflow in graphs, when we consider a graph as a simplicial 
complex of dimension 1. 
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2 PRELIMINARIES 



2 Preliminaries 



2.1 Flow in a network 

There are many equivalent ways to define the objects needed to state our problem. We will 
work with the following: 

Definition 2.1. A network is a pair (G, c) such that 

i) G — {V, E) is a. finite simple directed graph 

ii) V^V'U{s,t} 

in) for any v E V, {v, s) ^ E and (t, v) ^ E. 
iv) c:E — ^ K+ 

We call s and t the source and the sink respectively. Condition iii) means that there are no 
edges into the source, and no edges out of the sink. 

Definition 2.2. For any simple directed graph we define the incidence Junction 
(p:V X E — > {-1, 0, 1} as follows: 

{1 : e = {v,x) 
-1 :e = {x,v) 
:else 

Definition 2.3. A flow on a network [G, c) is a function / : E — > M_)_ such that: 

i) for any edge e, /(e) < c(e) 

ii) for any vertex v ^ s,t conservation of flow holds: 

J]/(e)0Ke)=O 
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Definition 2.4. For a flow / on a network {G, c), we define the net flow. 

I/I :=5^/(e)0(.,e) 

I /I is the total amount flowing out of the source. 



2.2 The problem of MAXFLOW 



Given a network {G, c) the MAXFLOW problem is to find a flow / of maximum net flow. 
Theorem 2.5. For any network (G, c) there exists a flow / of maximum net flow. 



Proof. Let m : [\E\] E he an enumeration of the edges of G. Let J-" be the set of feasible 
flows on (G, c). The map ip : ^ Ac M'^I, / ^ [f{m{i))]i is a bijection between J-' and 
a subset of M'^L Let J-'* be the image of J-' under this map. From definition 2^, the edge 
capacity constraints imply that J-"* is bound, and the flow conservation constraints imply 
that it is closed, hence J-'* is compact. The map f*cJ^*^ \'^~^{f*)\ € M is a linear map, 
hence continuous. As it is defined on a compact set, it achieves a maximum value, say at 
/max- ^/'"^(/max) is then a flow of maximum value. 

□ 



The previous theorem shows that, in fact, MAXFLOW is a linear programming problem, the 
most important results of which can be proved with LP theory. We discuss this formulation 
in detail in what follows. 

Definition 2.6. Let G = (V, E) be a simple directed graph, its incidence function and 
(G, c) a network. Let n := \V\, m := \E\, v : [n] ^ V he an enumeration of the vertices 
and e : [m] E he sua enumeration of the edges. We define the incidence matrix $^^ewith 
respect to the enumerations v,e as 

^^vAhj) = Av{i),e{j)) 
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Table 1: Incidence matrix of Gi 



From now on, we suppose a network {G, c) has fixed enumerations v,e, of vertices and edges. 
We take n as the number of vertices, m as the number of edges and suppose that v{l) = s, 
v{n) = t, then we refer simply to the incidence matrix as $. 

Lemma 2.7. Given a network (G,c), the problem of MAXFLOW is equivalent to the fol- 
lowing LP problem: 

max $1 ■ X : = : ImX < c* : x > 

where $ij is the first row vector of the matrix $, $* is the matrix that results from $ by 
deleting the first and last rows, is the identity matrix of size m and c* := [c(e(z))]j is the 
vector of edge capacities. 



Proof. Follows from the definition 2.3 and the proof of theorem 2.5 



□ 

Definition 2.8. For a linear program (called the primal problem) 

maxfc^x) ; Ax < b ; x > 
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the dual program is defined as 

min(6"^l/) ; A^y > c ; y > 
y 

We will compute the dual program of MAXFLOW. 

Definition 2.9. Let (G, c) be a network. A cut is a partition of V into two disjoint subsets 
(5', 5') such that s G 5 and t E S. Let i5' := {e G : e = (s, s) s e 5* , s e 5}. We define 
the capacity of the cut (5*, S) 

C{S,S) = J2<^) 

eeE' 

C{S, S) is the sum of the capacities of the edges directed from S to S. We say an edge e 
traverses the cut ii e — {s, s) s & S s & S 

Lemma 2.10. For any cut {S,S') 

\f\^f{S,S')-f{S',S) 
where f{X, Y) is the sum of the values of / at the edges directed from X to Y. 

Proof. 



I/I = E /(^)'^(^' ^) = E E /(^)'^(^' ^) = E E /(^)'^(^' ^) 
- E E/(^)'^(^'^)+ E E/(^)'^(^'^)+ E E/(^)'^(^'^) 

V ✓ V ^ V ' 



=f(S,S') -f{S',S) 



□ 



Lemma 2.11. Given a network (G, c), for any fiow / and any cut (5", S) 

\f\<C{S,S) 

The problem of MINCUT is to find a cut of (G, c) of minimum capacity. Two of the most 
important results are the following 



8 



2.2 The problem of MAXFLOW 



2 PRELIMINARIES 



Theorem 2.12. MAXFLOW=MINCUT. This means the net flow of a maximal flow is equal 
to the capacity of a minimal cut. 

In order to prove these results, we will see that the dual program of MAXFLOW is a 
relaxation of the MINCUT problem and use the following: 

Theorem 2.13. (weak duality). Let x* and y* be feasible solutions to a primal problem and 
its dual, respectively, then 

c X < y 
Proof. For yi, 2/2 > define the function 



fl-d/i, 1/2) := maxc x + y^{b- Ax) +y2X 

X 

Clearly, for any feasible x* and yi,y2 > 0, c^x* < g{yi,y2)- 



Vlb -.c^-yjA + y^^O 
00 : else 



Rearranging terms we have 

g{yi: Vt) = max(c'^ - y^A + yl)x + y^h 

X 

where yi G and 2/2 G I^"- Then we have 

9{yi,y2) = 

Now minimizing over yi, ^2 > we have, for any feasible x* 

c^x* < min g{yi,y2) 

2/1 ,2/2 >0 

By the previous observation this is equivalent to min^^^ y2>o ulb subject to — yiA + y2 = 0. 
We see that — yjA + y'^ — ^ — yjA — —y'^ < and this is equivalent to a single 
variable y > such that — y'^A < =^ < y'^A and this is the dual problem, as we 
wanted to show. 

□ 
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Theorem 2.14. (strong duality). If the primal problem has an optimal solution x*, then 
the dual problem also has an optimal solution y* and 

c X —by 

To find the dual of our problem, we state it in standard form 



; x>0 

Where 0„ is the zero column vector of length n. Then we find the dual to be 





$* 




On-2 


X ; 


— $* 


X < 


On-2 








C 



(1) 



mm 

j/gK2(n-2)+m 



On-2 
On-2 
C 



-I T 



After further inspection this is equivalent to n — 2 unrestricted in sign variables, one for each 
vertex v{i) ^ s,t called Vi, and m variables, one for each edge Cj such that 



min Je ; [ Im][v,ef>^T. ; e>0 

[i),e]eIR"-2+'" 

These restrictions translate to the following set of inequalities: 

Vi-Vj + ek>Q: v{i) ^ s, t v{j) ^ s, t e{k) = {v{i), v{j)) : i < j 
-Vi + ek>l : v{i) ^ s,t e{k) = (s, v{i)) 
Vi + Ck > : v{i) s,t e{k) — {v{i),t) 
efc > VA; e {l,...,m} 

We can define vi — —1 and Vn — and we write all the equations in the form 

Vi-Vj + ek>0: e{k) = {v{i),v{j)) : i < j 



(2) 



10 



2.2 The problem of MAXFLOW 



2 PRELIMINARIES 



Lemma 2.15. For any cut {S, S) there exists a feasible solution of ([2j) such that the value 
of the function at this feasible solution equals the capacity of the cut. 

Proof. Let {S, S) be a cut. Define = 1 if and only if e{k) = (s, s) with s E S, s E S and 
Cfe = in any other case, f « = — 1 if and only if v{i) & S and f j = in any other case. Then 
c^e = C{S, S) and it is straightforward to check that the restrictions hold. 

□ 

As a corollary we get 



Corollary 2.16. Lemma (2.11) 



Lemma 2.17. For an optimal solution [t'*,e*] of ^ there exists a cut {S,S) such that 
C{S,S) < Je*. 

Proof. Let x ^ ["IjO] be a random variable with uniform distribution. Define a random 
variables for each edge by 



eik) 



X 



1 ■.v*<x< V* 
: else 



Note that this assignment defines a random cut. If v* < v* then P(t>* < x < "^j) = 
min{l, V* — V*} < V* — v* < el then by the restrictions of the problem, we get 

E{C{S,S)) = 5^CfeE(e(A;)J < J^c^e^ = cV 

As the expected value of the random cut capacity is less or equal to the optimal value of the 
problem, there exists a cut of capacity less or equal to the optimal value. 

□ 

This proves that the dual of MAXFLOW is in fact a relaxation of MINCUT and we get, by 
strong duality 



Corollary 2.18. Theorem (2.12) 
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3 MAXFLOW algorithms 



3.1 The Ford-Fulkerson algorithm 

L.R, Ford Jr. and D.R. Fulkerson devised a polynomial time algorithm to compute a maxi- 
mal flow first published in 1962 |lj. We introduce some new concepts needed to describe the 
algorithm, and prove some general facts about it. 



Definition 3.1. Given a network {G, c) and a flow / we define c, f : V x V ^ 



c{u, v) 



c{u, v) : {u, v) G E 
: else 



f{.u,v) 



f{u,v) 







{u,v) e E 
{v,u) e E 

else 



Lemma 3.2. For a given flow / on a network {G, c), for any fixed vertex u ^ s,t 
Remark 3.3. In fact, it is equivalent to define a flow as a function / mapping V x V to 



the reals such that the equation in definition 3.2 holds, and such that f{u,v) = —f{v,u) for 
every pair of vertices. From now on we refer to a flow in this sense, and we refer to /, c as 
/ and c whenever it does not cause confusion. Under this new definitions we have that the 
capacity of a cut can be written as 

C{S,S')= J2 ^("'^) 

u£S,v£S' 



and lemma (2.10) translates to 



1/1= E 



for any cut {S, S') 
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Definition 3.4. For a network {G, c) and a flow / we define the residual capacity of a pair 
of vertices {u, v) as c/(m, v) = c{u, v) — f{u, v). Any sucli pair witli residual capacity greater 
than zero is called a residual edge. Note that the residual capacity is always greater or equal 
to zero. We define the residual graph Gf as the graph with vertex set that of V{G) and edge 
set the set of residual edges. 

Definition 3.5. Given a flow / on a network (G, c), an augmenting path is a directed path 
on Gf from source to sink. 

Lemma 3.6. a flow / is maximal if and only if there is no augmenting path on Gf. 



Proof. Suppose there is no augmenting path on Gj. Let S be the set of vertices v such that 
there exists a directed path from s to w in Gf. Let S' = V\S. {S,S') is then a cut. By 
definition of S, we have that c/(m, v) = ^ c{u, v) = f{u, v) for any u E S, v E S' then we 
have, following remark (3.3) 



C{S, S') 



I/I 



'xes,v£S' 



So / is a maximal flow by Theorem (2.12) or Theorem(2.13). 



Now suppose there is an augmenting path {vq = s,Vi, . . . , Vk-i, Vk = t). Let A = min{c/(f j, fj+i) : 
z = 0, . . . , /c - 1} > 0. Define F : V x V R as F{vi, Vi+i) = f{vi, fj+i) + A, F{vi+i, Vi) = 
/(f j+i, fj) — A and F{u, v) = f{u, v) on any other pair of vertices. One can easily check that 
F is a flow, and that \F\ = \ f\ + A so f is not a maximal flow. 

□ 

Now this theorem is the basic result needed to state the Ford-Fulkerson algorithm. Starting 
with the zero flow, as long as there exists an augmenting path with respect to such flow, we 
can increase the value of the flow by A as defined in the above proof. 

Lemma 3.7. If {G,c) is such that c{u,v) G Z then the algorithm terminates. 

Proof. At each step of the algorithm, the value is increased by A > 1 so a maximal flow is 
reached after a finite number of steps. 



□ 
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As corollaries we get 

Corollary 3.8. If the capacities of a network are integers, then the value of the maximal 
flow is an integer and there exists a maximal flow with f{u, f ) G Z for every edge {u, v). 

Corollary 3.9. If the capacities of a network are rational numbers, then the algorithm 
terminates. 

In fact there are examples of networks with irrational capacities such that the algorithm never 
terminates, moreover, the value of the flow in each step does not converge to the actual value 
of the maximal flow, so our algorithm must have as a condition that the capacity is at least 
a rational valued function. Then, the running time of the algorithm depends on the way the 
augmenting paths are chosen. There are many ways to find an augmenting path, like the 
shortest augmenting path or the largest bottleneck (value of A) augmenting path, that lead 
to a polynomial time algorithm. 



3.2 The Goldberg-Tarjan algorithm 



The Goldberg-Tarjan algorithm |2] is another polynomial time algorithm with a different 
approach to the problem of finding a maximal flow. Instead of increasing the flow along 
augmenting paths, it starts with a preflow, which is a function on V x V which satisfies 
excess of flow at each vertex, and then pushes excess flow to edges closer to the sink. Next 
we formalize these concepts following Goldberg- Tarjan's article |2,|. 

Definition 3.10. Given a network {G, c) a preflow is a function f : V x V satisfying: 

i) fiu,v) < c{u,v) 

ii) f{u,v) = -f{v,u) 

iii) for any vertex x s, J2v&v fi^^ v) >0 
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3.2 The Goldberg-Tarjan algorithm 3 MAXFLOW ALGORITHMS 

Definition 3.11. for a network (G, c) and given a prefiow / on the network, we redefine 
the residual capacity of {u,v) G V x V as Cf{u,v) = c{u,v) — f{u,v). If Cf{u,v) > we call 
such pair a residual edge. We define the residual graph as the directed graph having vertex 
set V and edge set the set of residual edges. 



Note there are similarities with definition (3.4) but in this definition we are working with a 
prefiow rather than a flow. 

Definition 3.12. the excess flow at a vertex x & V is defined as Yliv<^v /(^''^) — 0- 

Definition 3.13. given a a valid labeling on a network (G, c) is a function d -.V Z>oU{oo} 
such that d{s) = n, d(t) = and d{v) < d{w) + 1 for every residual edge (m, f). 

It can be shown that for any vertex f , if d{v) < n then d{v) is a lower bound on the distance 
from f to t in the residual graph and if d{v) > n then d{v) — n is a lower bound on the 
distance to s in the residual graph [2]. This labeling of the vertices permits the algorithm 
to push excess flow to vertices that are estimated to be closer to the sink and, if needed, to 
return flow to vertices estimated to be closer to the source. 

Definition 3.14. a vertex v is called active if d{v) < oo and e{v) > 0. 

Now we define the basic operations, push and relabel, that the main algorithm uses. 

Push. Let {v,w) be such that v is an active vertex, Cf{u,v) > and d{v) = d{w) + 1. 
Define 5 = min(e(t'), Cf{v, w)) > 0. Redefine f{v, w) := f{v, w) + 6, f{w, v) := f{w, v) — 6, 
e{v) := e(f ) — 6 and e{w) := e{w) + 6. 



Relabel. Let v be an active vertex such that for any w eV, Cf{v,w) > ^ d{v) < d{w). 
Redefine d{v) := mm{d{w) + 1 : (f , w) residual edge}. 

As initial prefiow we take the function / such that for any v G V, f{s, v) = c{s, v) and zero 
everywhere else. It is readily checked that this is a prefiow. As an initial labeling of the 
vertices we take d{s) = n and zero everywhere else. As long as there is an active vertex v, 
either an operation of push or relabel is applicable to v. When there are no more active 
vertices the algorithm terminates, and the prefiow becomes a fiow, and in fact, it is maximal. 
Details of the proof of correctness and termination of the algorithm can be found in [2] . We 
show only correctness assuming termination. 
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Lemma 3.15. If / is a preflow and d is any valid labeling for / then the sink t is not 
reachable from s in Gj. 

Proof. Suppose (s = Vq, . . . ,Vk = t) is a path from s to t in the residual graph. Clearly 
k < n. Now is a residual edge for every i. So by definition of valid labeling 

d{vi) < d{vi+i + 1) so we have d{s) < d{vi) + 1 < . . . < d{t) + k = k this contradicts the fact 
that d{s) = n. 

□ 



Now recall lemma (3.6) 



Theorem 3.16. If the algorithm terminates and c/ is a valid labeling for / with finite labels, 
then / is a maximal flow. 

Proof. At the end of the algorithm there are no active vertices, as the labels are finite, it 



means that all vertices have zero excess, so / is a flow. By lemma (3.15) and lemma (3.6) 
this flow is in fact maximal. 

□ 



One important remark about this algorithm is the fact that it always works (it terminates 
and it is correct) no matter what type of capacity function we are dealing with. The Ford- 
Fulkerson fails to terminate in some cases where the capacity function is not rational. It 
is also important to note that the algorithm relies only on local operations, that means 
the operations depend and modify only parameters related to a small part of the graph, 
this allows a parallel implementation of the algorithm that takes advantage of multicore 
processors. A special implementation of such algorithm terminates after Oin^m) steps. 



3.3 Hochbaum's pseudoflow 

Dorit Hochbaum's pseudoflow algorithm f3] is an algorithm with a different approach to 
the maximum flow problem. Instead of directly finding a maximum flow, it first solves the 
maximum blocking cut problem, then a maximum flow is recovered. Although the most 
complicated of the three, it is also the most efficient. We follow [3): 



16 



3. 3 Hochbaum 's pseudoQow 



3 MAXFLOW ALGORITHMS 



Definition 3.17. A pseudoflow / on a given network (G, c) is a function f : V x V ^ 
such that 

i) f{v,w) = —f{w,v), V(f,w) eV xV 

ii) < f{v,w) < c{v,w), V(f,w) G E 

The concept of pseudoflow drops the conservation of flow constraint, preserves the capacity 
constraint on the edges of the graph and the antisymmetry constraint on V" x V". We define 
the residual capacity and residual graph in the same manner we did with flows and preflows. 



Definition 3.18. For a directed, weighted, simple graph G = (V, E) with weights for each 
V E V and arc capacities c{v, w) for each {v, w) E VxV,we will define Gst as a directed graph 
with vertex set Vst = VU{s,t},edge set Est = EU A{s)LI A{t) where A{s) = {(s,w) : > 0} 
and A(t) = {{v,t) : < 0} and arc capacities c{s,v) = w^, c{v,t) = —w^ and the other 
arc capacities left unchanged. Starting from Gst = (Kt, Est) we define the extended network 
G^^* as the graph obtained from Gst by identifying s, t as a single vertex r and adding the 
edges (r, t>), (t>, r) for every v G Kt\{s, t}. 



We define the excess of flow at a vertex e{v) as in definition ( 3.12). 



Now we consider a pseudoflow / on Ggt and a rooted spanning tree with root r, T of G"^^* 
such that 

i) / saturates all arcs in A{s) U A{t) 

ii) For every arc in E\T, f is either zero or saturates the arc. 

iii) In every branch all downward residual capacities are strictly positive. 

iv) the direct children of r are the only vertices that do not have zero excess. 

Definition 3.19. a spanning rooted tree with root r of G*^^* that satisfies the previous 
conditions is called a normalized tree. Note that this is an undirected graph. 

A child r, of r is classified as: 
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i) Strong if e(f ) > 

ii) Weak if e{v) < 

A vertex v is called weak or strong if it has a weak or strong ancestor, respectively. 

As we mentioned earlier, Hochbaum's algorithm solves first the maximum blocking cut prob- 
lem, which we state next: 



Problem: For a directed, weighted graph G = {V,E) with vertex weights w{v) for each 
vertex v, and arc capacity function c(a, b) defined for every (a, b) G E, find S C V such that 

surplus{S) = — c(a, b) 

bes 

is maximum. Such a set is called a maximum surplus set and {S, S) is called a maximum 
blocking cut. 



The key is to find the relation between a maximum blocking cut in G and a minimum cut 
in Gst- Given by the following lemma: 

Lemma 3.20. {s} U is the source set of a minimum cut in Gst if and only if {S, S) is a 
maximum blocking cut in G. 



This is proven in [3] following an article by Radzik |4J. The following lemma, also found on 
the article [s], is fundamental for the correctness of the algorithm. 

Lemma 3.21. For a normalized tree T and pseudofiow / on Gst saturating A{s) and A{t) 
and a set of strong vertices S, if the residual capacity of any edge (a, b) with a E S and b G S, 
Cf{a, 6) = is zero then S* is a maximum surplus set and {S, S) is a maximum blocking cut. 



For a normalized tree T if the set of strong vertices S satisfies the condition in lemma (3.21 ) 
the tree is called optimal 



The algorithm starts with a normalized tree related to a pseudofiow / on Gst- There are 
multiple choices of such a tree. We will start with a simple normalized tree. It corresponds 
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to a pseudoflow / saturating A{s) and A{t) on Gst- In this normalized tree every vertex 
in V forms an independent branch. The set of strong vertices are those adjacent to the source. 



By lemma (3.21 ), it is desirable to reduce the residual capacity from strong to weak vertices, 
therefore, with each iteration of the algorithm, a residual edge from S* to is chosen, this 
is called a merger arc(edge). If such an edge does not exist then the tree is optimal and the 
set of strong vertices form a maximum blocking cut. If there is one, then such edge becomes 
a new edge of the tree and the edge joining the root of the strong branch to r is removed 
from the tree. Then the excess of the root of the strong branch is pushed upwards until it 
reaches the root of the weak branch. Note that this path is unique. 



It is not always possible to push the total of the excess along an edge. If there is an edge, say 
(a, h) that does not have enough residual capacity to push the excess then such edge is re- 
moved (split) from the tree, a (the tail of the edge) becomes the root of a new strong branch 
with excess equal to the excess pushed minus the residual capacity of the edge. This is done 
in such a way so that the property that only roots of branches may have nonzero excess is 
maintained through the running of the algorithm. The remaining excess at h continues to 
be pushed in the same fashion until it reaches the root of the weak branch or until it reaches 
another edge that does not have enough residual capacity and the process is repeated. This 
process assures that the tree is normal at the end of each iteration. 



Termination of the algorithm follows from the next lemma: 



Lemma 3.22. At each iteration of the algorithm either the total excess of the strong vertices 
is strictly reduced or the number of weak vertices is reduced. 



Proof. Recall that from the properties of definition (3.19 ) we have that all downward residual 
capacities of edges are positive. After appending a merger edge to the tree and removing 
the edge joining the root of the strong branch to r, the path from to the weak branch 
becomes an upward path with positive residual capacity at each edge of the path, then some 
positive amount of excess arrives at the weak branch that is being merged. Then either 
some positive amount of excess arrives at the root of the weak branch and the total excess is 
strictly reduced, or there is some edge in the weak branch without enough residual capacity. 
In this case the edge is split and the tail of such edge becomes a strong vertex. Note that 
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if some weak vertex becomes strong in this fashion, then all of its children, including the 
former strong branch, becomes strong. Then if such operation takes place, the number of 
weak vertices is strictly reduced. 

□ 

Now let M"*" = C{{s}, V) be the sum of capacities in A{s) and M~ be the sum of capacities in 
A{t) then by the final comment in the previous lemma we see that any iteration that reduces 
the total excess is separated from another iteration of such type by at most n iterations. 
Then it follows immediately for integer capacities that 

Corollary 3.23. The complexity of the algorithm is 0{nM^) 

Now as the problem is symmetrical on s and t we find that by reversing all directions of 
the edges of the graph and interchanging s and t we get an equivalent problem so it follows 
again that for integer capacities 

Corollary 3.24. The complexity of the algorithm is 0{n * minlM"*", M~}) 



Correctness of the algorithm follows from lemma (3.21 ) as at the end of the algorithm there 
are no merger arcs left. 



Now in order to solve our initial problem we have to recover a maximum flow from the 
pseudofiow and maximum blocking cut obtained after the algorithm terminates, as it is 
not guaranteed that the pseudofiow becomes a flow after termination. In what follows we 
describe how to recover such maximum flow. 

Definition 3.25. An s — t path-flow on a network [G, c) is a flow / on (G, c) such that the 
edges carrying a strictly positive amount of flow form an s — t path on [G, c). A cycle-flow 
on (G, c) is a flow on (G, c) such that the edges carrying a strictly positive amount of flow 
form a directed cycle on (G, c). 

Theorem 3.26. (Flow decomposition) Let / be a flow on (G,c), then / can be decom- 
posed as the sum of at most m s — t path-flows and cycle-flows. 

Proof. Suppose / is such that |/| > then there is some G E such that f{s,vi) > 0. 

If f 1 = if: we have a directed s — t path and we define a flow /q carrying an amount of flow 
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f{s,vi) on such a path and zero everywhere else. If vi ^ t then there exists some edge 
(fi,f2) with some positive amount of flow as a result of conservation of flow. In this way 
we construct an s — t path (we may suppose it has no loops) and we define the flow /q as 
carrying an amount of flow equal to the minimum of the flow over the edges of this path and 
zero everywhere else, it is readily checked that this is a feasible flow. /' = / — /o is again a 
feasible flow where f{u, f ) > and f'{u, v) = 0. Using the same argument for /' we arrive at 
a flow /* = / — /o — • • • /fc with zero net flow. If /* is not the zero flow, then analogously to 
the previous argument we may construct a cycle on the graph and define a cycle-flow as the 
minimum over the flow of the edges on the cycle and zero everywhere else, this is a feasible 
flow fk+i and /* — fk+i has some new edge with zero flow. We continue in such fashion and 
arrive at f — f o — ... fh = so f = fo + ■■■ fh where fi is either a path-flow or a cycle-flow, 
/i < m as at least the flow on one edge becomes zero in each step. 

□ 

Lemma 3.27. (see js]) For any strictly strong node there exists a residual path either to 
the source or to some strictly weak node. 

In order to use the flow decomposition theorem first we have to consider a network (G', c') 
related to {G, c) such that the preflow / becomes a feasible flow. This is done by considering 
a super source s and supersink i, adding edges (s, s) U {(s, f ) : v is strictly weak} with flow 
and capacity equal to the deficits on such vertices, and edges (t, t) U {(f , t) : v is strong} 
with capacity and flow equal to the excesses on such vertices. The flow on any other edge 
has the same value as the preflow. This function is now a feasible flow on the network with 
source s and sink t. 



To get a feasible flow on the original network, we have to get rid of excesses at strong nodes 
and deficits at strictly weak nodes. For any strong vertex Vs, as long as /(fg, t) > we have 
that (t, f s) G Ef is part of the residual network. Hence by lemma (3.27) we have a residual 
path from t to s that contains the edge {i, Vg). Increasing the flow on such path by an amount 
of S equal to the minimum over the residual capacities of the path, actually decreases the 
excess of Vs by the same amount. After one such step, either the vertex Vs arrives at zero 
excess or this process can be repeated by lemma (3.27). This is a process analogous to flow 
decomposition on the reversed graph. After termination there are no vertices other than t 
with positive excess. 



In the same fashion, the remaining flow is decomposed until positive deficits at strictly weak 
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vertices are disposed. This must be done via t as it is the only vertex sending a positive 
amount of flow to t. After termination all vertices except s and t have nonzero excess. 
Deleting s, t from the graph leaves us with a feasible flow on Ggt- 

Corollary 3.28. A maximum flow can be recovered from an optimal normalized tree with 
pseudoflow f. 



Proof. For an optimal tree we have Cf{S, S) = that is, there are no residual edges directed 
from strong to weak vertices. Hence, following the previous argument, excesses at strong 
vertices can be disposed using only paths traversing strong nodes. Now there are no edges 
directed from a weak to a strong vertex with positive flow, as otherwise the reverse edge 
would have residual capacity greater than zero, a contradiction. So by the proof of theorem 



(3.26) the remaining deflcits at weak vertices are disposed using only paths traversing weak 



vertices. It then follows that Cf{S, S) = after recovering a flow / so that c(f , w) = f{v, w) 



for V G S,w G S and as consequence |/| = C{S,S). By lemma (3.20) and lemma (3.21) 
{S, S) is a minimum cut. This shows |/| is maximum. 

□ 
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4 Applications 

There are many not so obvious applications of maximum flow algorithms and results to 
different pure and applied topics, we show three interesting problems that can be solved 
using the previous results. 



4.1 Hall's Marriage Theorem 

Let G = (yuW,E) he a bipartite graph, where V nW = (J} and \V\ = \W\ = n. Label the 
vertices in \^ as . . . , f„, and the vertices 'm.W as Wi, . . . ,Wn- A perfect matching on G is 
a permutation a & Sn such that [f j, ifcr(i)] G E for every i = 1, . . . , n. 

Definition 4.1. Let S <Z V . N{S) := {w e W : [v,w] G E for some u G 5*} is the set of 
neighbors of S. 



We prove the following using the Maxfiow-Mincut theorem (2.12). 



Theorem 4.2. (Hall's Marriage Theorem) A perfect matching exists if and only if 

WS C V, \S\ < \N{S)\ 

Proof. Clearly such condition is necessary as a is injective. Suppose WS C V, \S\ < |A^(S')|. 
Now we construct an s — t network by directing all edges e G -E from V to W, adding a 
source s and sink t and appending the edges {(s, f ) : f G V} U {{w, t) : w E W}. We set the 
capacity of such new edges to 1, and the capacity of the original edges to n + 1. Let S" be a 
minimum cut on such network. We show that C{S, S) = n. C{S, S) < n as the cut S' = {s} 
has capacity n and C{S,S) is minimum. Now we show C{S,S) > n. Let X = S f\V . 
N{X) C W and if N{X) (t S then there would be an edge crossing the cut, of capacity n + 1 
so C{S, S) > n + 1 then N{X) C S' fl W. On the other hand, all edges traversing the cut 
are of the form (s, v') where v' G V\X or of the form [w', t) where w' E S Cl W. Then 

C{S,S)= J2 J2 c{w',t) = \V\S\ + \SnW\>n-\X\+N{X)>n-\X\ + \X\ 

v'ev\s w'esnw 
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So the capacity of a minimum cut is n. By theorem (2.12) and (3.8) there exists a maximum 
flow / of integer values and net flow |/| = n. As there are only n edges out of the source s 
and into the sink t, and they have capacity 1, they must be saturated. By conservation of 
flow and the fact that the flow is integer, for any v & V there exists only one w E W such 
that f{v,w) = 1. Again by conservation of flow and integrality, for any w E W there exists 
only one v E V such that f{v,w) = 1. This shows that the edges directed from V to W 
carrying a flow of 1 define a perfect matching on G. 

□ 

Corollary 4.3. There exists a polynomial time algorithm that finds a perfect matching on 
a bipartite graph. 



4.2 Counting disjoint chains in finite posets 

Definition 4.4. A finite poset P := {P', <) is a finite set P' together with a partial order 
< on P'. We say that P has or (1) if there exists an element x E P such that x < y or 
X > y for any y E P, respectively. A chain is a subset c := {xq, . . . ,Xn} C P such that for 
any two elements Xi,X2 either xi < X2 or X2 < Xi. 

Given a finite poset P, we say that a chain C is maximal if C U {x} is not a chain for any 
X E P\C. Clearly any maximal chain contains and 1. We say that y covers x in the poset 
if X < y and there exists no z such that x < z < y. We say that a set {Ci} of chains are 
cover-disjoint if whenever y covers x then {x, y} belongs to at most one chain Q. We would 
like to find a subset S of the set of maximal chains, such that 5* is cover-disjoint and such 
that 15*1 is maximum. 

One of the possible ways of doing this is to work in a greedy algorithm fashion, finding one 
of such chains and then repeating the process in the remaining part of the poset. We note 
that this may not lead to a partition of maximum size, as the example in Figure 2 suggests. 

Instead, we consider an associated network P^t where s = 0, t = l, V = P and (x, y) E E 
whenever y covers x . We define a capacity function with value 1 on every edge. 

Lemma 4.5. The maximum number of disjoint chains in P is equal to the net fiow of a 
maximum fiow in P^f. 
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Figure 2: Two sets of maximal chains in a poset 



Proof. We first show that given a set of disjoint chains Ci, . . . , we can find an associated 
feasible flow on Pst with flow value equal to A;. C^'jO, 1} U Cj forms a chain from to 1 and 
a directed path from source to sink in Pgt. Define a flow fi in Pgt as having value 1 on the 
edges of C'^ and zero everywhere else. This is a feasible flow. / = /j + . . . + is a function 
that satisfies conservation of flow and, as the chains were disjoint, it also satisfies capacity 
constraints so it is a feasible flow. Each fi saturates one edge leaving the source, hence 
I/I = k. This shows that k is always less than the value of a maximum flow. 



Now given a maximum flow / on Pst we construct a set of disjoint chains of size |/|. We may 
assume / has integer values by corollary (3.8). By theorem (3.26) we may write f = fi . . . fj 
where each /j is an s — t path. As the flow has integer values, so do /i, . . . fj. As the 
capacities are all equal to 1 each fi must have net flow equal to one and so the /j's do not 
intersect as they saturate all the path. Then {/i, . . . , fj} define a set of disjoint 0, 1 chains. 
Finally |/| = | A| + . . . + |/,| = j. 



□ 



4.3 Image segmentation 

The problem of segmenting a given image is that of defining a partition of the pixels as two 
sets, the foreground and the background, so that they form coherent regions. There are 
multiple other problems defined under the label of image segmentation. In the following we 
show how to define the problem and how to solve it using algorithms of flow optimization 
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in a network, following T.M. Murali's lecture notes |6j. Throughout this section we denote 
a directed edge as {v,w) and an undirected edge as [v,w]. 

We define a finite undirected graph G = (V, E) where V is the set of pixels of an image, 
V C X Z_(_ and the set of edges E comprises the set of neighbors for each pixel (x, y) G V. 
The set of neighbors of {x,y) is A''(x,y) := {{x + l,y), (x — l,y), {x,y + 1), {x,y — 1)} fl V. 
We define functions a : V [0, 1], f i-)- = probability that f is in the foreground, 
b : V ^ [0,1], V i-^ bj, = probability that v is in the background and a penalty function 
p : E ^ ]R+, [v,w] p[v,w] = penalty for defining v,w for defining v in the foreground and 
w in the background. 

Problem: Partition the set V as two sets A,B (foreground/background) such that the 
function 

s(v4,5) := + ^ 6^ - ^ piy^^i 

VGA weB [v,w]&E 

\An[v,w]\=l 

is maximized. The idea is that if > b^ it's preferable to set v as in the foreground and if a 
pixel V has most of it's neighbors defined as in the foreground, it is preferable to set v as in 
the foreground also. Such probabilities are given in the problem, however, different choices 
of such values may lead to better or worse results in the segmentation of the image. For 
instance if one is interested in isolating a small object in a big background, the best choice 
is to take higher values for the probability function a„. 

In order to construct such sets, one must define as foreground (background) vertices those 
with higher probability of belonging to the foreground (background), while reducing the total 
penalty of the boundary between foreground and background. We want to formulate this 
problem as a Mincut problem. In order to do this we have to overcome some difficulties, 
namely, that of working with an undirected graph rather than a capacitated network, and a 
function to be maximized rather than minimized. 

Lemma 4.6. Let Q = Ylvev ^" then 

s{A,B)=Q-^b^-^a^- p[^^^]=Q - s'{A,B) 

VGA wGB [v,w]gE 

\Anlv,w]\=l 
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where 

s'{A,B) = Y,bv + Y,^^+ Yl 

vi^A wGB [v,w]eE 

\An[v,w]\=i 

Then maximizing s{A, B) is the same as minimizing s'{A, B) 

Now we consider a directed graph G* = {V U {s,t}, E* U A{s) U A{t)] where E* = {{v, w) : 
G E}, A[s) = {{s,v) : V E V} and A(t) = {{v,t) : v G and we define a capacity 
function as c{s,v) = a^, c{v,t) = and c{v,w) = p[v,w]. We then have a network where 
the source (sink) is connected to each pixel with such edge with capacity equal to (K) and 
where each undirected edge [v,w] of neighbor pixels is replaced with two antiparallel edges 
(f , if), (w, v) both with capacity equal to p[f , w\. Then it follows immediately that for a cut 
{A,B) in such network we have 

C{A,B) = ^K + ^a^+ =s'(A,S) 

v(^A weB [v,w]gE 

\An[v,'w]\=i 



By lemma (4.6) arrive at next result: 



Corollary 4.7. A minimum cut {A, B) in G* solves the problem of image segmentation. 

There's only one difficulty left to overcome, as we must deal only with simple graphs, we 
must replace the set of antiparallel edges. We do this by adding for each pair of neighbor 
vertices v,w two new vertices Cyw,Cwy and replacing the antiparallel edges with the edges 
{v,Cyw),{cyw,w),{w,Cu,y),{cwv,v) sll with capacity equal to p[v,w]. It use readily checked 
that it is equivalent to find a maximum fiow on this new graph. 



Then we solve the problem by finding a maximum fiow on such network using either the 
Ford-Fulkerson or Goldberg- Tarj an algorithm and then recovering a minimum cut using 
theorem (3.6). 
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Figure 3: How to get a simple graph 

5 A Generalization of Maxflow 

We would like to define a more general optimization problem that reduces to the Maxflow 
problem on graphs and then try to generalize the optimization algorithms studied on previous 
sections. 

5.1 Preliminaries 

Definition 5.1. A simplex on a set 5" is a finite subset S' C S. 

Definition 5.2. A simplicial complex A on a set 5" is a set of simplices on S closed under 
taking subsets. Elements of A are called faces. Maximal faces (faces that are not subsets 
of any other face) of the complex are called facets. Elements of a simplex X are called its 
vertices. The dimension of a face X is defined as dim(X) = |X| — 1. The dimension of the 
complex dim(A) is defined as the maximum over the dimension of its faces. 

Definition 5.3. A simplicial complex A is called pure if all facets have the same dimension. 

Given a network (G, c) we can consider the graph G* that results from appending the edge 
(s, t) with infinite capacity and then, the problem of finding a maximum flow on the original 
network is equivalent to finding a maximum circulation on G*, that is a positive function on 
the edges of G* satisfying capacity constraints and flow constraints on every vertex. In this 
case the objective function |/| is the amount flowing through the edge with infinite capacity. 
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Definition 5.4. Let X — {vq, . . . , Vd} be a simplex with d> 1. Consider the set of orderings 
of vertices of X, 0{X) :— {[v^^o), ■ ■ ■ ;'i^o-(<i)] ■ e -S'd+i} modulo the relation [vq, . . . ,Vd] — 
[vo-(o), . . . ,t'cr(d)] ^ o" is even. This partitions the set in two equivalence classes that we call 
orientations of X. To choose an orientation for X is to choose one of such orientations, 
which we call the positive orientation and we say that X is oriented. We denote an oriented 
simplex as X = {vq, . . . ,Vd). 

Notation. For a simplicial complex X we let X^'^^ be the set of its d-dimensional simplices. 

Definition 5.5. For d > 1 Let Cd be the free abelian group over the orderings of elements 
in X^*^) modulo the relations [vq, . . . ,Vd] = [va{o), • • • , ^^(d)] ■v4> o" is even and [vq, . . . ,Vd] — 
~bo-(o)) ■ ■ ■ ) ''^aid)] ^ c" is odd. Co is defined in the same fashion but notice that the relations 
become trivial. 

Definition 5.6. The boundary operator 9^ : — > Cd-i is a homomorphism defined in the 
basis as 

d 

dd{X) = ^(-l)*[wo, . . . , Vi-i, Vi, Vi+i, ...,Vd] 

i=0 

where Vi means deleting such term. 

Elements of Cd are called d-chains. Elements of the subgroup ker{dd) are called ci-cycles. 
5.2 Higher Maxflow 

Definition 5.7. A d-dimensional network is a triple {X,T,c) where 

1. X is a simplicial complex of pure dimension d all of whose facets have chosen orienta- 
tions. 

is a distinguished oriented simplex of dimension d satisfying the source 

condition: 

• For every oriented d-simplex a which intersects T" in a (d — 1) -dimensional simplex 
T the signs of r in da and dT are opposite. 

3. c : X^'^) R+ is a function with c{T) = oo. 
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Remark 5.8. The source condition is a generalization to the assumption that on a network 
(G, c) every edge incident with the source is directed out of it and every edge incident with 
the sink is directed into it. The source simplex T is the generalization of an appended edge 
directed from sink to source with infinite capacity. 

Definition 5.9. A flow on a network (X, T, c) is a function / : X^'^^ — t- ]R_|_ satisfying the 
following properties: 

1. / is a weighted cycle, that is 

2. For every (i-dimensional simplex a we have < f{a) < c(cr) 

Remark 5.10. The condition that / is a weighted cycle is a generalization of the conserva- 
tion of flow condition. To see this we define the following: 

Definition 5.11. Let v = {vq, . . . ,Vd-i} E X^'^^^^ be a c/ — 1 dimensional simplex of a 
simplicial complex X of degree d. Fix an orientation [v] = [vq, . . . ,Vd\ of v. Let e G X^'^^ 



<P{v,e)= { -1 




1 '■ +[v] appears in de 
— [v] appears in de 
else 



This is, in fact, a generalization of the incidence function defined on section (2.1 ). Then the 
condition that / is a weighted cycle is equivalent to: 



for any d — 1 dimensional oriented simplex v 



J2 ^{v,e)f{e) = 

Definition 5.12. The amount carried by a flow / in a ((i)-dimensional network (X, T, c) is 
the number /(T). 

This comes from the fact that after appending the edge (t, s) to a network {G, c) with flow / 
one must define f(t,s) = \ f\ so that conservation of flow holds in all vertices including s,t. 



30 



5. 3 As an LP problem 



5 A GENERALIZATION OF MAXFLOW 



(HMax-Flow.) The higher max flow problem asks to find the maximum possible amount 
/(T) which can be carried by a fiow on a network (X, T, c). 

Remark 5.13. A (l)-dimensional network is a capacitated graph (T is the edge from t to 
s) and HMax-Flow reduces to Max-Flow on graphs. 

Remark 5.14. If all the capacities are 1 and X is a triangulated orientable (i-manifold and 
T is any top-dimensional simplex then every top-dimensional cycle is a fiow with /(T) = 1. 
This is one HMax-fiow. This follows directly from the definition of oriented manifold. 



5.3 As an LP problem 

Even in higher dimension, the problem can be stated as a set of linear equalities and inequal- 
ities in a finite dimensional vector space, so it can also be stated as a linear program as we 
did in section (2.2). We continue with the convention that \n\ = \X^'^~^^\ is the number of 
(d — 1) faces, and m = \X^'^'' \ — 1 is the number of facets without considering the source T. 
Then there are 2m fiow restrictions, two for each edge, and n fiow conservation restrictions, 
one for each vertex. As always we consider fixed enumerations of the {d — 1) faces Vi and of 
the facets Cj (suppose Vm+i = T). After fixing an orientation of the (d — 1) faces, we define 



the matrix [d]i 



(fj,ej), it has dimension n x m + 1. Let be the identity matrix of 



dimension m and J^, be such matrix with a zero column vector appended as the rightmost 
column. Then we can state the problem as : 



max Xm+i ; 



d 




' 0, " 




-d 


X < 


On 


; a; > 


Im^ 




C 





(3) 



where c G M™" is the vector q = cici] 



Computing the dual we find it can be stated as 



mm c e 



{?;,e)eM"+™ 

where f G M" and e G -R™. 



9 



T 



[v,e) > Cm+i, e > 
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5.4 Further examples and conjectures 

Example: Consider the three dimensional simplex. It is a simplicial complex on the set 
{1, 2, 3,4} where the elements of the complex are the faces of the simplex. 



4 




Figure 4: Oriented simplicial complex, the tetrahedron Sim2 

Take the orientation as the induced by the outward normal. With this orientation the 
oriented faces are [234], [124], [143], [132]. We choose the orientations of the one dimensional 
faces to be [12], [14], [24], [13], [34], [23]. The capacity on each facet is chosen to be 1. In this 
example we drop the commas, dealing only with integers less than 10. It is readily checked 
that [132] satisfies the source condition so we may take T = [132]. With the enumeration of 
facets and {d — 1) faces given by this order we have that the incidence matrix is given by: 






1 





-1 





-1 


1 





-1 


1 














-1 


1 


1 





-1 





1 








-1 



Table 2: Incidence matrix of Sim2 



32 



5.4 Further examples and conjectures 



5 A GENERALIZATION OF MAXFLOW 



Towards a good definition of a cut on a generafized c?-dimensional network we find feasible 
solutions of the dual of such problem. We are mostly interested in integer solutions to the 



problem. In this example we find both d and 







to be totally unimodular matrix so 



the existence of integer optimal solutions is assured. We suggest the following definition. 

Definition 5.15. A cut on a generalized (i-dimensional network X is a partition {S,S') of 
the d — 1 faces of X. 

Prom a cut we can construct a solution to the dual program of HMaxfiow. Let {S, S') be a 
cut. Assign to the first n dual variables A^, (where each one corresponds to a — 1 face v of 
the complex) the value if v e 5 and the value 1 if v e 5". For each facet a of the complex, 
define +(cr) be the set of d — 1 faces such that r induces their positive orientation and define 
— (o") in the same fashion. In the dual program, the variables that correspond to the d — 1 
faces are unrestricted in sign, while the variables that correspond to the facets rjfj must be 
nonnegative. There is one inequality in the dual program for each facet that can be written 
in the form 

ve+{<7) ve-{<7) 

for any a ^ T, and 

Aj, - ^ A^, + ?7t > 1 

v&+{T) fe-(T) 

So, given the cut and the assigned values to the d — 1 faces variables we can define the dual 
variables corresponding to the facets rj^ as 



for a and 
Vt = 



"(J^vG+ia) ■^v — J2jj&-(a) ■^v) : else 

-(E^e+(a) - E^e-(a) ^ - 1) : else 

In such way it is readily checked the solution is dual feasible. We define the capacity of a 
cut (5", S') as the value of the dual objective function at this solution, which is the weighted 
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sum of the capacities of facets a ^ T, with weights rjcr. In such way the capacity of a cut 
is always an upper bound to the value of a maximum flow, by weak duality. It is natural 
to ask whether the minimum of the capacities over all cuts equals the value of a maximum 
flow. We attempt to show this using an analogous probability argument as in section (2). 

As in section ^ we can extend the capacity function defined on X^'^^ to the set of all 
orientations of elements of X*^'^^ as c(x) = c(x) if x G X^'^^ and c(x) = if -X G X^*^) and 
f{x) = f{x) if X G X('^) and f{x) = -f{-x) if -x G X^"^). 

Definition 5.16. The residual complex Xj is the (multi) simplicial complex whose facets 
are those x G X^"^) U — X'^'^^ such that the residual capacity of x Cf{x) := c{x) — f{x) > 0. 
By definition, if the capacity function is not identically zero, the residual complex is a pure 
multisimplicial complex of dimension d. 



Recall from that from lemma (3.6), a flow is not maximal if there exists a simple (no loops) 
path from s to t in the residual graph. After appending the edge (t, s) this reduces to: a 
flow is not maximal if there is no simple cycle containing (t, s) on the residual graph. This 
motivates the following definition: 

Definition 5.17. Given a (X, T, c) a d dimensional network and a flow / on X, an aug- 
menting cycle is d-cjcle a such that cr = ^ . CjXj with Cj G Z+ and Xj G Xf and such that 
Xj = T for some i. 

Lemma 5.18. Let (X, T, c) be a network and /, /' be two feasible flows in such network. 
Then / + /' is positive and it is a weighted cycle. 



Proof. 

The fact that / + /' is positive is clear. □ 
Lemma 5.19. A flow is not maximal if there exists an augmenting cycle. 



Proof. Let A = ^ CiXi be an augmenting cycle with Xj G Xf and Cj positive. Let m = 
min{^c/(Xj)}. Define a flow as /(X) = /(X)+Cjm if X = Xj for some i, /(X) = f{X)—Cim 

if X = -Xi for some i and /(X) = /(X) else. Now define B = {X e X^'^) : X = X^orX = 
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J2 Kx)d{x)= Yl f{x)d{x) + Y,Kx)d{x) 



j\:ex('')\B 



= /(^)^(^) + 5Z ±Cim9(X) = md\y^c,XA=Q 

xexM xeB \ i J 

Capacity constraints hold by definition of m. As Xj = T for some i the value of the flow is 
strictly increased. 

□ 



We would like to prove the converse of this lemma to devise a first algorithm for Higher 
Maxflow optimization. 



Conjecture 5.20. The converse of lemma (5.19) holds in general 




Figure 5: Oriented simplicial complex, the double tetrahedron Sim2,2 



Example Consider two tetrahedra oriented by the outward normal, with a common facet 
T, over the vertex set {1, 2, 3, 4, 5} We find the incidence matrix to be 
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This matrix is also totally unimodular hence the problem has integer optima. In various cases 
we find this optima to be the sum of the minimum capacities over each tetrahedron. This 
optimum is achieved after finding two augmenting cycles corresponding to each tetrahedron. 
When the fiow is maximum there is no augmenting cycle on the residual path. 

Question 5.21. When is the incidence matrix of a d dimensional network totally unimod- 
ular? 



We will cite some important theorems that will give us some insight about the status of 



the question (5.21), and ultimately show that it is in fact false. We find certain family of 



networks where it holds. 



The map dd has a unique matrix representation with respect of chosen basis. We denote 
such matrix as [dd\- In fact it coincides with the incidence matrix we have defined. The 
kernel of dd is called the group of p-cycles and is denoted by Zd{X). The image of the 
map dd+i, denoted by Bd{X) is a subgroup of Cd called the (i-boundaries. We have that 
for any d dddd+i = under composition so that Bd C Zd- The rf-dimensional homology 
group is defined as the quotient Hd := Zd/Bd- For a simplicial couplex X we have that 
Zd{X), Bd{X), Hd{X) are finitely generated abelian groups and Hd{X) depends only on the 
homotopy type of X. 



For a subcomplex Xq C X we define the group of relative chains of X modulo Xq as the 
quotient Cd{X)/Cd{XQ) := Cd(X, Xq). The map dd induces a map 0^^'^°^ : Crf(X, Xq) — )■ 
Cd-i{X,Xo) which also satisfies df'^°^d^Jlf'^ = So we define the relative homology groups 
as Hd{X,Xo) = fcer(9f '^"V^^lCf"^) ■= Xo)/S,(X, Xq). 
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The following is a result by Dey-Hirani-Krishnamoorthy that gives a partial result to our 
question: 

Theorem 5.22. [5j For a finite simplicial complex triangulating a + 1 dimensional com- 
pact orientable manifold, [dd+i] is totally unimodular irrespective of the orientations of the 
simplices. 



The answer to question (5.21) is "not always". Consider the following counterexample: 



Counterexample: [s] For certain simplicial complex triangulating the projective plane, the 
matrix [82] is not totally unimodular. 

This might not be exactly a counterexample for our conjecture as we need to prove existence 
of some facet that may work as source. However if we consider the two dimensional sphere 
positively oriented, and a triangulation of such manifold with consistent choice of orienta- 
tions, then we may define any facet of such complex as the source, and then joining this 
complex to the triangulation of the projective plane by a vertex we find a submatrix of [^2] 
that is not totally unimodular. 



The following is a theorem also due to Dey-Hirani-Krishnamoorthy that characterizes totally 
unimodular matrices arising from a boundary operator. 

Theorem 5.23. |5j [9^ + 1] is totally unimodular if and only if Hj^{Y, Yq) is torsion free, for 
all pure subcomplexes Y, Yq of X of dimensions d + 1 and d respectively, where Yq C Y. 



The following theorem yields another family of simplicial complexes where total unimodu- 
larity holds, namely, the family of d dimensional complexes embeddable in M.'^. 

Theorem 5.24. |5| Let i^' be a finite simplicial complex embedded in ]R"'+^, then Hd{L, Lq) 
is torsion free for all pure subcomplexes Lq and L of dimensions d and d + 1 respectively. 



It would be helpful to find more general families of simplicial complexes where unimodularity 
holds. Using theorem (5.23) we may give an alternative proof of the fact that for graphs, 
the incidence matrix is totally unimodular. 

Theorem 5.25. For a directed graph G, the incidence matrix [di] is totally unimodular. 
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Proof. Suppose G is connected. Let Lq C Lhe two subcomplexes of G of dimension and 1 
respectively. We want to show that Hq{L,Lq) is totally unimodular. Now {L,Lq) is a good 
pair, this implies that Hq{L,Lq) ~ Hq{L/Lq) which is a torsion free group isomorphic to 
Z'^ where k is the number of connected components of L/Lq. The result follows by theorem 
( [K23| ). 

□ 

Definition 5.26. Let X be a simplicial complex. A pair of facets (F, F') of X is a leaf if for 
every face if of X we have F n H (Z F' . A simplicial tree is a simplicial connected complex 
X such that every subset of facets of X contains at least one leaf. 

Lemma 5.27. Let X be a pure simplicial complex of dimension 2, and suppose that for any 
pure sub complex of dimension 2 L C X, Hi{L) = 0, then [82] is totally unimodular. 



Proof. Let Lq C L be pure sub complexes of dimension 1 and 2 respectively. We have a long 
exact sequence in homology 

...^Hi{L)^H,{L,Lo)^Ho{Lo) 

By exactness of the sequence we have that Hi{L,Lq) — > Hq{Lq) is an invective map. As 
Hq{Lq) is torsion free so is Hi{L,Lq). The result follows by theorem (5.23). 

□ 



Lemma 5.28. Let X be a pure simplicial complex of dimension 2 such that there exist 
disjoint facets {Ti, . . . , T^} such that X\{Ti, . . . , T^} is a simplicial tree. Then for any pure 
subcomplex L C X of dimension 2, Hi{L) = 0. 



Proof. Let L be a pure simplicial complex of dimension 2 with facets {Fi, . . . ,Fi}. Define 
V as the pure simphcial complex with facets {Fi, . . . , F/}\{Ti, . . . , T}. By assumption, V is 
a two dimensional subcomplex of a two dimensional simplicial tree so it is also a simplicial 
tree. Any simplicial tree is contractible and homology groups are homotopy invariant. Now 
for any F facet of L that is not in V we have two cases: Every one dimensional face of F is 
in V or not. In the first case, after contracting V, F forms a sphere. In the second case F 
has at least a one dimensional facet not in V and F can be contracted to FnV. Hence after 
contracting V we see that L is homotopy equivalent to a wedge of spheres so that Hi{L) = 0. 

□ 
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Corollary 5.29. Let X be a pure simplicial complex of dimension 2 such that there exist 
disjoint facets {Ti, . . . , T^} such that X\{Ti, . . . , T^} is a simplicial tree. Then [^2] is totally 
unimodular. 



Conjecture 5.20 remains without answer. 
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